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Abstract 

We obtain sharp asymptotic estimates for hitting probabilities of a critical branch¬ 
ing Brownian motion in one dimension with killing at 0. We also obtain sharp asymp¬ 
totic formulas for fhe tail probabilities of the number of particles killed at 0. In the 
special case of double-or-nothing branching, we give exact formulas for both the hit¬ 
ting probabilities, in terms of elliptic functions, and the distribution of the number of 
killed particles. 


1 Introduction 

Branching Brownian motion is a stochastic particle system in which each individual particle 
moves along a Brownian trajectory, and at a random, exponentially distributed time inde¬ 
pendent of its motion is replaced by a random collection of identical offspring particles. 
The motions, gestation times, and offspring numbers of different particles are condition¬ 
ally independent, given the times and locations of their births. Thus, conditional on the 
event that at time t there are Zt particles at locations xi,X 2 , - ■ ■, xzt, the law of the post- 
t evolution is identical to that of Zt mutually independent branching Brownian motions 
started by individual particles at the locations x*. A formal construction of the process is 
outlined in section 2 below. 

The process Zt that records the total number of particles at time f is a continuous-time 
Galton-Watson process: see [AN72], ch. 2 for the basic theory of these. Branching Brownian 
motion is said to be supercritical, critical or subcritical according as the mean of the offspring 
distribution is greater, equal or less than 1. In the critical and subcritical cases, the particle 
population eventually dies out, with probability one, provided the population starts with 
only finitely many particles; in the supercritical case, however, there is positive probability 
that the population blows up, that is, Z* —)■ cx) as t ^ oo. Thus, the questions that are 
germane to the supercritical case are different from those of interest in the critical case. 

It has been known since the work of McKean [M75] that supercritical branching Brow¬ 
nian motion is intimately related to the behavior of solutions to the Fisher-KPP equation. In 
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particular, this equation governs the cumulative distribution function u{t, x) = F^{Mt < x) 
of the position Mt of rightmost particle at time t. Using this fact, McKean gave a proba¬ 
bilistic proof of the Kolmogorov, Petrovsky, and Piscounov [KPP37] theorem, which asserts 
that the solution of the KPP equation with Heaviside initial data stabilizes as a traveling 
wave of velocity ^/2. Subsequently, Bramson [B78] used the connection with supercritical 
branching Brownian motion to obtain sharp estimates for the center of the wave, and Lal- 
ley and Sellke [LS87] showed that the limiting traveling wave w{x) can be represented as 
a as a mixture of extreme-value distributions. 

When the branching mechanism is critical or sub-critical, a more natural object of study 
is the random variable 

M = sup Mt, 
t>o 

the rightmost location ever reached by a particle of the branching Brownian motion. Criti¬ 
cal branching Brownian motion has been proposed as a model for the spatial displacement 
of alleles without selective advantage or disadvantage, and in this context the distribu¬ 
tion of M plays an important role (see, for example, [CG76], [S76], [SF79] and references 
therein). Sawyer and Fleischman [SF79] proved that if the offspring distribution has mean 
1, positive variance and finite third moment, then the tail of the distribution of M satis¬ 
fies the power law 

g 

P(M > x) ~ „ „ as X —)• oo. (1.1) 

cr^x"' 

Modifications of branching Brownian motion and branching random walk in which the 
laws of reproduction and/or particle motion depend on particle location arise in various 
contexts. See, for instance, Lalley and Sellke [LS88], [LS89] and Berestycki et al. [BBS13], 
in which particle reproduction is allowed only in certain favored regions of space; Kesten 
[K78], Aldous [A], Addario-Berry and Broutin [ABll], Aidekon, Hu, and Zindy [AHZ13] 
and Maillard [M13] where particles are killed upon entering the half-line (—oo,0], and 
Berestycki et al. [BBHM15]; and Lalley and Sellke [LS92] and Korostelev and Korosteleva 
[KK03], [K04], [KK04], where particles move according to spatially-inhomogeneous diffu¬ 
sion laws. In the articles [ABll], [AHZ13], [M13], and [BBHM15], the branching law is su¬ 
percritical, but particle production is balanced by the killing in (— 00 ,0] so that M„/n —)• 0. 

This paper will focus on the modification of critical branching Brownian motion (that 
is, where the mean number of offspring at reproduction events is 1) in which particles 
are killed upon reaching the interval (— 00 ,0]. Clearly, the number Zt of particles alive at 
time t in this process is dominated by the corresponding random variable for the critical 
branching Brownian motion with no killing, and so = 0 eventually, with probability 1. 
Furthermore, the distribution of the maximal particle location M is dominated by that of 
the maximal particle location in critical branching Brownian motion with no killing, and 
so the results of Sawyer and Fleischman [SF79] imply that for any e > 0 and any initial 
particle location y > 0, 

P^M > x} < 

for all sufficiently large x. 
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It is by no means evident, however, that the tail behavior should be the same as for 
branching Brownian motion with no killing. In fact we will prove that when the branching 
process is initiated by a single particle at a location y > 0 near zero, the tail follows a power 
law with exponent 3 rather than 2. In particular, we will prove in Theorem 6.1 that for each 
fixed y > 0, 

P^(M > x) ~ as X —>• oo. 

where Cs > 0 is a constant depending on the offspring distribution but not on x or y. On 
the other hand, we will show that for initial particle locations y = sx whose distances from 
the killing zone are proportional to the target x, the exponent of the power law reverts to 
2; in particular, there exists a continuous function 04 ( 5 ) of s G (0,1) such that 

P*^(M > x) ~ as X ^ 00. 

Furthermore, we will show that in the Moranian case, where the offspring law is double-or- 
nothing, the tail probability (M > x) can be explicitly written as a Weierstrass P—function. 

All of these results will be deduced from an analysis of a boundary value problem in the 
variable y satisfied by the hitting probability P^(M > x). 

Also of interest is the total number N of particles killed at 0. For supercritical branch¬ 
ing Brownian motion with particle drift and killing at 0, Maillard [M13] and Berestycki 
et al. [BBHM15] have, under various hypotheses concerning the drift and the reproduc¬ 
tion mechanism, obtained sharp estimates for the tail of the distribution of N. For critical 
branching Brownian motion with killing, T. Y. Lee [L90-1] proved a conditional limit theo¬ 
rem for the distribution of N given that A > 1: in particular, he showed that as the position 
y of the initial particle —)• 00, so that P^(A > 1) —)• 0, the P^— conditional distribution of 
N/py{N > 1) converges to a non-degenerate limit distribution. (See also [L90-2] for a 
time-dependent analogue.) 

We will study the distribution of N for critical branching Brownian motion with killing 
at 0 under a fixed P^. In section 7 we will show that, for offspring distributions with mean 
1, positive finite variance cr^, and finite third moment, 

k 

Y.jFy{N>j)^CyVk as k ^ 00. 
i=i 


Under certain additional hypotheses on the offspring distribution, we will show that the 
distribution of N obeys a power law with exponent 3/2, thus verifying a conjecture of 
Professor Jian Ding, and in addition, we will show that the N obeys an asymptotic local 
limit theorem. In particular, we will prove that 


P2^(A > A:) ~ 


Cry 

A;3/2’ 


and 

~ 

where Cr,C% > 0 are constants depending on the offspring distribution. Finally, in the 
Moranian case, we will give in Theorem 7.9 an explicit formula for the tail distribution of 

N. 
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2 Construction and Monotonicity Properties 


Branching Brownian motions with initial particle locations at points y G M+ can be con¬ 
structed on any probability space that supports countably many (i) independent standard 
Wiener processes Vh*; (ii) independent, identically distributed unit exponential random 
variables T*; and (iii) independent, identically distributed random variables Lj all dis¬ 
tributed according to the prescribed offspring distribution. We dub this construction the 
discrete Brownian snake, as it is the natural discrete analogue of Le Gall's Brownian snake: 
see [L99] for details. 

The construction proceeds by using the random variables {Ti}j>o to construct a Galton- 
Watson tree. This construction is standard: see [AN72]. If the offspring distribution has 
mean 1, as we shall assume throughout, then the resulting Galton-Watson tree is almost 
surely a finite, rooted tree with vertices arranged in generations, beginning with the root 
at generation 0. To each vertex v is attached one of the random variables Li, with Lq 
attached to the root; for each vertex v the random variable Li determines the number of 
offspring vertices. The random variables Li can be attached to vertices in any number of 
different ways, the most common being the breadth-first rule, in which the values Lj are 
read successively from the stack generation-by-generation, left-to-right. 

Given the realization of the Galton-Watson tree, we attach unit exponential random 
variables Tj and standard Wiener processes W* to the edges of the tree in such a way that 
the index i = i{v) matches the index of the random variable Li attached to the lower vertex 
V of the edge (the incident vertex with higher generation number). The random variable 
Ti attached to an edge determines the real time elapsed between reproduction events, and 
the Wiener process W* determines the displacement of the particle in real time from its 
position at the last reproduction event. Thus, the particles alive at (real) time t are in one- 
to-one correspondence with the vertices v of the tree such that 

Ti(u,) < t ^ Tii^w)] 

w<v w<v 

here the symbols < and < indicate the ordering of vertices w along the geodesic path in 
the tree from the root to v. The spatial position of the particle represented by vertex v at 
time t is 

y + Y, + w^t - Y T^M)■ 

w<v w<v 

Observe that these rules yield a simultaneous construction of branching Brownian motions 
from all initial positions y. It is evident from this construction that the distribution of the 
maximum position M attained by a point of the branching Brownian motion is stochasti¬ 
cally monotone in the initial position y. 

Branching Brownian motion with killing at 0, or more generally with killing at any 
point z < 0, can be constructed using the same marked tree as for branching Brownian 
motion with no killing. The rule is simple: once a trajectory along an edge enters (—oo, z], 
the tree is pruned at that point. This leaves a subtree of the original Galton-Watson tree 
in which certain edges (those corresponding to particles that are killed at 0) are cut. The 
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vertices of this subtree represent particles of the branching Brownian motion with killing 
at z. Thus, the set of particles alive in the branching Brownian motion with killing is a 
subset of the set of particles in the branching Brownian motion with no killing, which we 
will henceforth refer to as the enveloping branching Brownian motion. 

This construction makes it obvious that the distribution of M is dominated by that for 
branching Brownian motion with no killing at z, and that if Z 2 < zi then the distribution 
of M for branching Brownian motion with killing at zi is stochastically dominated by that 
for branching Brownian motion with killing at Z 2 - Furthermore, the implied inequalities 
among the cumulative distribution functions are strict: for instance, if W 2 {x) and wi{x) 
are the tail distributions of M for branching Brownian motions with killing at Z 2 < zi < 0, 
respectively, when both are initiated by a single particle at 0 (that is, Wi (x) is the probability 
that M > x) then 

< W2{x). {1.1) 

To see this, observe that there is positive probability that a branch will be pruned when 
there is killing at zi but not when the killing is at Z 2 , and that this branch will extend 
in such a way that it gives rise to a particle that reaches location x. Finally, branching 
Brownian motions with killing at z converge as z —)■ — oo to branching Brownian motion 
with no killing. Thus, for any x > 0, 

lim Wz{x) = Woo{x), {1.1) 

Z^OQ 

where Wz{x) is fhe probability that M > x for branching Brownian motion with killing at 
—z and Woo{x) is the corresponding probability for branching Brownian motion with no 
killing (both with initial particles located at 0). 

It should be obvious that minor variations of the construction just outlined can be used 
to build a variety of related processes. One that will prove useful in certain of the argu¬ 
ments to follow is branching Brownian motion with freezing, in which particles that reach a 
target point 0 (or, more generally, a closed set B) are frozen in place, ceasing all motion and 
reproduction thereafter, but not dying. In a critical branching Brownian motion with freez¬ 
ing of particles at location 0, eventually all existing particles will be frozen at 0; moreover, 
the number Nt of particles frozen at time t is the same as the number of particles killed at 
0 up to time t in the corresponding branching Brownian motion with killing at 0. 

Henceforth, we shall assume that all branching Brownian motions are critical and that 
the offspring distribution has positive, finite variance and we shall denote by 

OO 

^{z) = Y,^iL = k)z^ (2.3) 

k=0 

the probability generating function of the offspring distribution. 


3 Product Martingales and Differential Equations 

The key to our analysis will be the fact that hitting probabilities and related expectations 
for critical branching Brownian motion, viewed as functions of the initial point y, are gov- 
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erned by a nonlinear second-order differential equation. This differential equation is well 
known, but since we will have occasion to consider expectations of complex-valued ran¬ 
dom variables, we shall spell out the boundary value problems in detail. 

Say that a sequence / : C is multiplicative if it is a geometric sequence of the form 

/(n) = z"' for some z £ C. For any A G (0, oo], let be the law of a branching Brownian 
motion with initial point y G [0, in which particles are frozen upon reaching either 0 or 
A. For i = 0 and i = A define 

Ni = number of particles frozen at i. (3.1) 

Both A^o arid Na are almost surely finite, since only finitely many particles are born in the 
course of a critical branching Brownian motion. Clearly, Na = 0 when ^ = 00 . 

Proposition 3.1. if /, <7 : Z+ —)• C are bounded, multiplicative sequences, then the function 
(f{y) = W [f{No)g{NA)] satisfies the second order differential equation 

= (p{y) - T'((p(y)) for all y G (0,^). (3.2) 

In the special case where d. = 00 and /(n) = So{n) this was stated and proved by 
Sawyer and Fleischman [SF79], and this proof was subsequently cited by Lee [L90-1]. But 
the proof in [SF79] seems to have a gap: the derivation of the differential equation relies 
on the smoothness of the function ip{y), but to prove this the authors quote the version 
of WeyTs Lemma given in [M69] to conclude that a weak solution must be C°°. We do 
not understand this argument, as WeyTs Lemma, in the form stated in [M69], applies only 
to linear parabolic differential operators, while the differential operators in [SF79], section 
2, and in our Proposition 3.1 are nonlinear. Therefore, we will sketch another approach 
to the proof of Proposition 3.1 that uses an interesting class of product martingales. (Sim¬ 
ilar martingales for supercritical branching Brownian motion were used in [LS 88 ] and 
[N 88 ]). Let h : [0, —)• C be a function bounded in absolute value by 1, and denote 
hy Xi{t),..., Xz{t){t) the locations of the particles alive at time t (including those frozen at 
one of the endpoints 0, A) in a branching Brownian motion with freezing at 0 and A; define 

z(t) 

Y{t) = Yh{t) = Wh{Xi{t)). (3.3) 

i=l 

Proposition 3.2. If h{y) satisfies the differential equation h” = h — T'(/i) in the interval (0,d.) 
then Y{t) is a bounded martingale, relative to the standard filtration for the branching Brownian 
motion, under ¥y,for any y G [0, A]. 

Proof of Proposition 3.1 (Sketch). Given Proposition 3.2, we proceed as follows. Fix /, g, and 
let /i : [0, j 4] —)• C be the unique solution to the boundary value problem 

lh" = h-^(hy, 

h{0) = /(i), 

h{A) = g{l). 


6 


(When A = oo, the boundary condition should be replaced by h(A) = g(0) = 1.) The exis¬ 
tence and uniqueness of solutions follows by standard arguments in the theory of ordinary 
differential equations; we omit the details.^ By Proposition 3.2, the process Y(t) defined 
by (3.3) is a bounded martingale, and so for any t < oo, 

EyY{t) = y(0) = hiy). 

But for all sufficiently large t, all particles will be frozen at either 0 or A, so eventually Y (t) 
coincides with f{NQ)g{NA). (For this the product structure of the martingale is essential.) 
Therefore, by the bounded convergence theorem, 

h{y)=EyfiNo)giNA). 


□ 

Proof of Proposition 3.2 (Sketch). By the Markov property, it suffices to show that for any 
initial configuration of particles y = {yi,y 2 ,... ,ym) the expectation E^y (t) is constant in 
time. Since each of the m particles engenders its own independent branching Brownian 
motion, the expectation E^y (t) factors as 

m 

Eyy(t) = ]jE^W(t); 

i=l 

consequently, it suffices to prove that for any y G [0, A] the expectation E^^y (t) is constant 
in time, and for this it is enough to show that 

But for this another conditioning shows that it is enough to prove that the derivative is zero 
at f = 0. This can be accomplished by a routine argument, by partitioning the expectation 
into the expectations on the events that the initial particle reproduces or not by time t 
and using the fact that h is bounded and (7^ and satisfies the differential equation h"/2 = 
/i-T'(/i). □ 

4 Weierstrass' V Functions 

In the special case of double-or-nothing branching (the Moranian case), the probability 
generating function of the offspring distribution is the quadratic function 'I'(s) = |(1 + 
s^). In this case the differential equation (3.2) reduces, as we will show, to the differential 
equation of the Weierstrass V—function. For a given period lattice 

L = {moj + riLo, m, n G N} , 

^At any rate the argument is routine in the case where /(I) and gr)!) take values in the unit interval [0,1]; 
in this case existence follows by a routine phase-portrait analysis for the associated first-order system, using 
the nonnegativity of the forcing term'I>(/i). When/(I) andg(l) are complex-valued, however, other methods 
must be used. See the proof of Lemma 7.14 in section 7 below for a proof in the case needed for the theorems 
on the distribution of the number of killed particles. 
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where oj and Cj are nonzero complex numbers whose ratio is not real, Weierstrass' V func¬ 
tion with period lattice C is the meromorphic function on C defined by 



See [K84] or [MM99] for expositions of the basic theory. Clearly, (4.1) defines a doubly- 
periodic function of z whose periods are the elements of fhe lattice C. It is also evident 
from (4.1) that P—functions with proportional period lattices are related by a scaling law: 
in particular, for any /3 7 ^ 0 and any laffice C, 

V^cil^z) = ^Vciz) for all z G C. (4.2) 

It is known (cf. [K84] or [MM99]) that the restrictions of Vc and its derivative to 
a fundamental parallelogram are branched covers of the Riemann sphere C of degrees 2 
and 3, respectively, and so for all but three exceptional values re G C the equation Vc{z) = 
w has two solutions zi,Z 2 in each fundamental parallelogram, and V'^{zi) = —V'^{z 2 )- 
Furthermore, the function Vc{z) satisfies the differential equation 

V'ciz)^ = AVcizf - g 2 {C)Vc{z) - gsiC), (4.3) 

where the constants g 2 {C) and gs (£) are given by the Eisenstein series 

g2iC) = QO ^ 

Sr{C) = 140 ^ i 

For any two complex numbers A, B such that — 27 / 0, there exists (cf. Proposition 
III. 13 in [K84]) a lattice C such that 


<72 (iC) = A and (4.4) 

gsiJO) = B. 

Proposition 4.1. Let A and B be two constants such that A^ — 27B‘^ / 0, and Let u{z) be a 
function on an open interval J c M with derivative u'{x) f Ofor all x G J that satisfies the 
differential equation 

ufzf = Mzf - Au{z) - B. (4.5) 

Then for some lattice L and some a G C, 

u{x) =Vc{x + a) for all x G J. (4.6) 

Proof Without loss of generality, assume that 0 G J and that u'(0) f 0. The differential 
equation (4.5) implies that in some neighborhood of x = 0, for one of the two branches of 
the square root function, 

u'{x) = ■\/4tt(x)3 — Au{x) — B (4.7) 
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Since u'(0) / 0, the right side of this equation is a Lipshitz continuous function of u(x) 
for X near 0, and so the Picard-Lindelof theorem guarantees that the equation (4.7) has a 
unique solution with initial value m( 0 ). 

Let £ be a lattice such that equations (4.4) hold. Because the Weierstrass P—function is 
a double covering of C, there exist two arguments a, a' G C such that Vc{oi) = Vc{a') = 
ri(0), and for one of these (say a) it must be the case that = m'( 0). Since the functions 

Vc{x + a) and Vc{x + a') both satisfy the differential equation (4.5), one of them (say 
Vc{x + a)) must also satisfy (4.7). By the Picard-Lindelof theorem, the equation (4.6) must 
hold in J. 

□ 

The connection between the differential equation (3.2) and the Weierstrass P—function 
is easily explained. If h{z) = z^, then the forcing term in (3.2) is quadratic, and so after a 
rescaling (3.2) can be written in the form 


u''{y) = 6u{yf. 

Multiplying both sides by u'{y) and integrating yields 

iu'iy)f=4u{yf + C, 


(4.8) 


(4.9) 


where C is a constant of integration. This is the characteristic equation for a P—function 
whose period lattice satisfies fl' 2 (£) = 0 . 

Proposition 4.2. The Weierstrass function u = Vc satisfies the differential equation (4.9)/or some 
C G C \ {0} if and only if the period lattice is of the form 

£ = {moj -|- (4.10) 

for some uj / 0; furthermore, C > Oin (4.9) if and only if the lattice has the form (4.10) with 

(4.11) 


a; = Icale"/®. 


In this case, u has real poles at integer multiples o/\/3|a;|, and takes only real values on R; fur¬ 
thermore, its only zeros in the fundamental parallelogram are at -\/3|a;|/3 and 2-\/3|a;|/3, and u is 
strictly increasing on (2\/3|t<;|/3, ■\/3|a;|). 

Proof. The Eisenstein series for the lattice £ = {muj + noj} can be written as 

* 

( 72 (£) = 60 ^(mw -|- nuj)~^ = 60 a;“^G 4 (^) and 

m,n 
* 

gsiC) = 140 y~^(mcc -|- nu))~^ = 140(4;“®G6(O 


(4.12) 


where ^ = w/tc is the ratio of two fundamental periods and the sum is over all pairs of 
integers except (0,0). By convention, the periods are ordered so that > 0; with this 
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convention, G 4 and Gq are modular forms of weights 4 and 6 (cf. [K84], section III.2). By 
the residue theorem for modular forms (cf. [K84], Proposition 111.2.8), any nonzero mod¬ 
ular form of weight 4 has precisely two zeros in the closure of the standard fundamental 
polygon of the modular group, at the points and Therefore, any 

Weierstrass function u = Vc that satisfies the differential equation (4.9) must have period 
lattice of the form (4.10) (as the choices and lead to the same lattice). 

The lattice (4.10) is invariant under rotation by vr/S (that is, C = and so by 

averaging over the six rotations one finds that 

G6(e"®/3) ^ ^ > q. (4.13) 

m^n 


Consequently, if gsiC) = —C < 0 then £ must be of the form (4.10) for some u such that 
cj® < 0 , that is, w is a positive multiple of a primitive 12 th root of unity. Thus, in the case 
gsiC) = —C < 0 the lattice £ must have the form (4.10) with lo = 


Assume now that £ is of the form (4.10) for some oj satisfying (4.11). Then by the 
addition law for the elliptic curve = 4x^ — 53 (cf. [K84], section 1.7; see especially Problem 
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Vc(V3\uj\/3) = rc{2V3\co\/3) = 0. 


Since Vc has degree 2 , it has only two zeros in a fundamental parallelogram, and by equa¬ 
tion (4.9) the derivatives V'^{V3\co\/3) = —T’£(2-\/3|t<;|/3) must be the two square roots of 
C. It is easily seen that the unique solution of (4.9) with initial conditions u{yo) = 0 and 
u'{yo) > 0 must be strictly increasing, with increasing derivative, on any interval {yo,yi) 
on which the solution u is well-defined and finite. This implies that V'^{V3\oj\/3) is neg¬ 
ative, and hence P(.(2\/3|ca|/3) is positive. It then follows that u is strictly increasing in 
(2V3|a;|/3,V3|.;|). 

□ 


Remark 4.3. The case where C = —g^{C) = — 1 in equation (4.9) is known as the equian- 
harmonic case; cf. [AS72] for further information. In the equianharmonic case the period 
lattice is of the form (4.10), but with to > 0, i.e., the lattice is of fhe same form as in the 
case where (7 = 1 but rotated by —tt/O. Call the case where C = —gz{C) = -|-1 the anti- 
anharmonic case; then by the scaling law, the P—functions for the equianharmonic and the 
anti-anharmonic cases are related by 

V z) = e~'^'-/^VEAH{z) for all z G C. 

Thus, mapping properties and special values of the P-function in the anti-anharmonic 
case can be read off from those for the equianharmonic case, which have been extensively 
tabulated. 

As far as we know, the occurrence of the P—function in critical branching processes 
was first observed by the first author in [L09], sec. 1.8. However, [L09] mistakenly asserts 
that the differential equation (4.9) with (7 > 0 falls into the equianharmonic case, and 
consequently the formulas in [L09], sec. 1.8 are off by factors of and 
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5 Distribution of M: Moranian Case 


In this section we consider the Moranian case [SF79], where the number of offspring is 
either zero or two, each with probability In this case the probability generating function 
is 'l'(t) = ^ + For 0 < y < X define 

w,(y) = P^{M > x} (5.1) 


to be the probability that the maximum position M attained by a particle of the branching 
Brownian motion initiated by a particle at y, with freezing of particles at 0, will exceed x. 
The function (fxiy) = 1 — Ux{y) is of the form covered by Proposition 3.1, so it satisfies the 
differential equation (3.2) with T'(z) = (1 + 2;^)/2, and consequently Ux satisfies 


<{V) = 2 


\ + Ux{.y)f - (1 - Ux{y)) 


= Ux{yf 


(5.2) 


Theorem 5.1. For branching Brownian motion with Moranian offspring distribution and killing 
at 0, the tail distribution function Ux{y) = P^(M > x) is given by 


Ux{y) = QVc^iy + 2ujx/2>), 


(5.3) 


where Vc^ (z) is the Weierstrass V function with period lattice 


Cx 




m—=e 

V3 


a/3 


(5.4) 


for some Ux > 0. The positive period ojx is uniquely determined by the boundary condition 


6T’/:.(a^ + 2a;a;/3) = 1. 


(5.5) 


Remark 5.2. The value of ujx can be computed numerically, by exploiting the fact that the 
inverse of the Weierstrass P—function w = Vc{z) is given by the elliptic integral 


z = 





dt 

92{F)t - gsiC) 


(5.6) 


For lattices of the form (5.4), we have g 2 {F) = 0. Consequently, the boundary condition 
(5.5) implies that 



dt 

\/4t3 - g^{£) 


(5.7) 


This equation (5.7) determines gz{F-), and hence, using the identities (4.12), the value of Wx- 
For X = 1, the values are 


g^{Ci) = -0.023786 • • • and wi = 9.88285 • • • . 

The large X dependence of Wx on x will be further clarified below, in Corollary 5.4. 
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The proof of Theorem 5.1 will rely on the uniqueness theorem for solutions of the differ¬ 
ential equation for the P—function (Proposition 4.1). For this, it will be necessary to know 
that u'^{y) / 0 for any y G [0, x\. 

Lemma 5.3. The function y u^iy) has strictly positive derivative u'^{y) on the interval [0, x]. 

Proof. Since Ux{y) is a cumulative distribution function, it is non-decreasing, and so its 
derivative must be nonnegative. Moreover, the differential equation u"{y) = Uxiy)"^ im¬ 
plies that the derivative is increasing at every y where Ux{y) > 0. It is easily seen that 
a branching Brownian motion started at any y > 0 has positive probability of putting a 
particle at x, so Ux{y) > 0 for all y G (0, x]. 

It remains to show that ri^(O) > 0. The differential equation u" = v? can be rewritten 
as the autonomous system 

u' = V, 
v' = uf 

The vector field in this system is clearly Lipshitz continuous, so solutions to the initial 
value problem are unique. Since u = r; = 0 is the unique solution with initial conditions 
u{0) = u(0) = 0, it follows that any non-constant solution of u" = vf satisfying tt(0) = 0 
cannot have derivative ri'(O) = 0. □ 

Proof of Theorem 5.1. Set Uxiv) = \ux{y)', then (5.2) becomes 

Ux{y)" = &Ux{yf, 

the differential equation encountered earlier in (4.8). The integrated form is (4.9). By 
Lemma 5.3, the derivative u'x{y) is strictly positive on y G [0,x], and so the same is ob¬ 
viously true of Ux{y). Hence, by Proposition 4.1, Ux must coincide with a translate of a 
P—function, and so for each x > 0 there exists a unique period lattice and a unique G C 
such that 

Ux{y) =Vc{y + ax) for all y G [0, x]. 

Since there is no linear term in the equation (4.8), the period lattice C must be of the form 
(5.4). 

Finally, the boundary condition rtx(O) = 0 and (5.3) imply that is a zero for Vc{z). 
Since Ux{y) is increasing in y and Vc{z) is doubly periodic, the constant Ux must be the 
larger zero of Vc{z) in (0, ojx), and in particular, by Proposition 4.2, 

Ux = 2a;a;/3. (5.8) 

□ 

All of the P—functions that occur in Theorem 5.6 are scaled versions of Vcx ■ By (4.2), if 
the positive periods wi and Ux of the lattices Ci and Cxr respectively, are related by 

Ax , (5.9) 
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then 


rcAz) = xlVcA^^z). 


(5.10) 


It is obvious that lim 3 ;^oo = 00 , because the function Ux{y) is an increasing function 
on (0, x), and hence cannot have a positive period smaller than x. Although the depen¬ 
dence of ujx on X is not linear, it is asymptotically linear, as the next corollary shows. 

Corollary 5.4. 

lim —= lini = 3. (5.11) 

x-foo X x^oo X^x 

Proof. The scaling law (4.2) and the boundary conditions for the functions Ux and ui imply 
that 

VcA‘^uji/ 3 + Xxx) = (5.12) 

Since (6A^)“^ —)• 00 as x —)• 00 , it follows that 2a;i/3 -|- XxX converges to wi, as this is the 
smallest positive pole of Vci ■ The result now follows from the equation (5.9). □ 

Remark 5.5. By exploiting the fact that Vci {oJi — z) ^ 1/as z ^ Q, one can obtain from 
the equation (5.12) the sharper approximation 

cOx = 3x + 3/\/6 -|- 0 ( 1 ) as X ^ 00 . (5.13) 

The scaling laws (5.10) and the period asymptotics (5.11) now combine to provide the 
large—X asymptotic behavior of the hitting probability function Ux{y). 

Theorem 5.6. For branching Brownian motion with Moranian offspring distribution and killing 
at 0, the tail distribution function Ux{y) = P^(M > x) of the maximum attained position M 
satisfies 

lim x^Uxiy) = Cl ■ y for each y > 0, (5.14) 

x^cyD 

where Ci = 6cfVjr^{2uji/3) = 33.0822 • • • and ci = cai/3 = 3.29428 • • • are constants that do 
not depend on x or y. Furthermore, for each fixed 0 < s < 1, 

lim x^Ux{sx) = C' 2 (s) (5.15) 

X^OO 

where C 2 {s) = 6cfP/;j(2a;i/3 -|- sci). 

Proof of Theorem 5.6. By equations (5.3) and (5.2), the function Vci satisfies the second- 
order differential equation 

'P'Uz) = 6VcAzf. (5.16) 

Using the abbreviation ax = 2ojxl3 and the fact that ^^^(ai) = 0 (cf. Proposition 4.2), it 
follows by taking successive derivatives that 
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Lemma 5.3 implies that V^^{ai) > 0, since this is proportional to tti(O). Consequently, 
since —)• 0 as x —oo, Taylor expansion around the point oi yields 

Ux{y) = GVc, {ax+ y) = Q^^l'Pci (ai + Ky) 

= 6A2 (v'c,{ ai)Xxy + (ai)A^/ + • • • ) 

= QV'c,{a^)\ly + 0 {\% 

Therefore, by Corollary 5.4, 

lim x^Uxiy) = 6 clVc^{ai)y = Ciy. 


and (5.14) follows. 

To prove (5.15), notice that for y = sx, we have 

x^Ux{y) = Qx'^Vc^ (ax + y) = Q{Xxx)‘^Vci (ui + s(A3;x)) 
By the continuity ofVci{z) on the interval (oi, cci), it follows that 

lim x‘^Ux{y) = Qc\Vci{ai + sci) = (3*2(s). 


□ 


6 Distribution of M: General Case 

In this section we will show that the asymptotic formulas (5.14) and (5.15) extend to branch¬ 
ing Brownian motions with arbitrary mean 1 offspring distributions with finite third mo¬ 
ments. The main result is as follows. 

Theorem 6.1. Assume that the offspring distribution has mean 1, positive variance cr^, and finite 
third moment. Then for each y > 0, the probability Ux{y) = P^{M > x} that a particle of the 
branching Brownian motion with initial particle at location y reaches location x satisfies 

lim x^Uxiy) = C 3 ■ y, (6.1) 

X^OO 

where €3 = Cij . In addition, for each fixed 0 < s < 1, 

lim x^Ux{sx) = 674, (6.2) 

X^OO 

where C 4 = Here Ci, 672 are the constants in Theorem 5.6. 

This theorem will be deduced from Theorem 5.6 by comparison arguments for differen¬ 
tial equations. The strategy is similar to that used by Lee [L90-1]: the key is that for small 
values of Ux{y), the forcing term h{ux{y)) in the differential equation 

u" = h{ux) (6.3) 
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(which follows from Proposition 3.1 as in the Moranian case) is well-approximated by the 
quadratic function a‘^Ux{yY- To see this, let T'(z) be the probability generating function 
of the offspring distribution. If the offspring distribution has finite third moment, then by 
Taylor expansion 

T'(l - z) = T'(l) - T''(l)z + O(z^) as z ^ 0. 

By hypothesis, the offspring distribution has mean 1 and positive variance 0 < < oo, so 

T''(l) = 1 and T'"(l) = Consequently, as z —)• 0, 

h{z) = 2 [^(1 -z)-{l- z)] = a^z^ + 0{z^). (6.4) 


Our arguments will use the following comparison principle for solutions to differential 
equations and inequalities. This is a minor modification of the Comparison Lemma in 
[L90-1]; because the result is standard and its proof involves only elementary calculus, we 
shall omit it. 


Lemma 6.2 (Comparison Principle). Let vi,V 2 be positive functions on an interval [yi, 2 / 2 ] such 
that for some constant a > 0, 

Vi{y) — avi{y)‘^ < 0 and (6.5) 

v'iiy) - av2{yf > 0 

for all yi <y < 2 / 2 . If 

vi{yi) > V 2 {yi) and (6.6) 

vi{y 2 ) > ^’ 2 ( 2 / 2 ), (6.7) 


then 


vi{y) > V 2 {y) for all 2/1 < 2/ < 2 / 2 - 


( 6 . 8 ) 


Next, we record several monotonicity properties of the functions Ux{y) = P^{M > x}. 

Proposition 6.3. The function Ux{y) is strictly decreasing in x and strictly increasing in y, and 
tf^(O) is strictly decreasing in x. 


Proof The monotonicity of Ux (y) in y and in x follow directly from the construction of the 
branching Brownian motion outlined in section 2. 

To prove that tt^(O) is strictly decreasing in x, recall that for each x > 0 the function 
Ux satisfies the differential equation u" = h{u), together with the boundary conditions 
^ia:(0) = 0 and Ux{x) = 1. Hence, by the uniqueness theorem for differential equations, if 
0 < xi < X 2 then f because otherwise the functions Uxfy) would be equal 

for all y € [0, xf, which is impossible because Ux^ (a^i) < 1- 

Thus, to complete the proof it suffices to show that if 0 < xi < X 2 then u'x^ (0) < 
u'x^ (0) is impossible. But since Ux-^ (0) = Ux 2 (0) = 0, if u'x^ (0) < u'x^ (0) then for all y 
in some neighborhood (0,e) we would have Ux^{y) < Ux 2 {y)- This would contradict the 
monotonicity of Ux{y) in x. □ 
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To study the behavior of Ux (y) near y = x, we introduce the function 


Wx{t) = Ux{x - t) (6.9) 

Observe that Wx{t) is the probability that a branching Brownian motion with killing at 
— {x — t) and initial particle at 0 will produce a particle that reaches location t. The con¬ 
struction in section 2 shows that Wx (t) is strictly monotone in both x and t, and that 

lim Wx{t) = Woo{t) (6.10) 

X—>-CxD 

where Woo{t) is the probability that a branching Brownian motion with no killing started 
at location 0 will produce a particle that reaches location t (cf. equation (2.2)). The conver¬ 
gence (6.10) holds uniformly for t in any finite interval [0,i*]. The function Woo{t) is the 
same as the function p{x) in the d = 1 case studied in [SF79], who proved that 

6 1 

rcoo(i) = + O(^) asf-s^oo. (6.11) 


To prove Theorem 6.1, we must determine the behavior of Wx{t) for large t, and in 
particular for t within distance 0(1) of x. The basic strategy will be as follows. For any 
e > 0 there exists < (X) so large that Woo{t*) < e. This implies that Wx{t:^) < e, or 
equivalently Ux{x — t^) < e, for all large x. Thus, in the interval [0, x — t*] the function 
Ux will be bounded above by e, and so h{ux) will be well-approximated by the quadratic 
function a'^Ux- The analysis of sections 4-5 shows that the differential equation (3.2) with 
h{u) = Cu^ admits an exact solution in terms of a Weierstrass P—function, so it will follow 
from the comparison principle above that in the interval [0, x — f*] the function Ux will be 
trapped between two such P—functions. By taking e —)■ 0, we will obtain sharp asymptotic 
approximations to Ux- 

Taylor expansion oih{z) shows that for all 0 < 5 < 1 there exists e = e((5) > 0 such that 
if 0 < Ua; (y) < e, then 


0-2(1 - 5)ul{y) < h {ux{y)) < o-2(l 5)ul{y). (6.12) 

On the other hand, (6.11) and the mono tonicity of rcoo(i) imply that for each e E (0,1), there 
exists ts such that Woo{te) = e and Woo{t) < e for all t > t^. By (6.10), for all 0 < y < x — 4, 

Ux{y) = Wx{x - y) < Woo{x - y). 

Therefore, (6.12) applies for all 0 < y < x — fg. Define 

r]{x,e) = Ux{x - te)] 

then for each e > 0 the function r]{x, e) is increasing in x, and so 

ri{x,£) =Wx{ts)^ Wco{ts) = e asx^oo. (6.13) 
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Corollary 6.4 (Pinching). Let te and r]{x, e) be as above, and define 


a+ = iT^(l + J) > 0 
a_ = a‘^{l — 5) > 0. 


If Uxiy),uf (y) and (y) satisfy the boundary value problems 


then 


' u'xiy) = h{uxiy)) 

Ux{0) = 0 

_ Ux{x - te) = r]{x,s) 


' uf'iy) = cLiU^'^iy) 

u±(0) = 0 

, utix - te) = y{x,e), 


ufiy) < uxiy) < {y) 


for all 0 < y < X — te- 


(6.14) 


(6.15) 


(6.16) 


Proof. This is an immediate consequence of the comparison principle (Lemma 6.2). □ 

The differential equations (6.15) for the functions uf are, except for the constants a±, 
identical to the differential equation (4.8) for the P—function. Consequently, they are re¬ 
lated by a simple scaling law. 

Lemma 6.5. Suppose that the functions uf{y) satisfy (6.15), and set 

utiv) = ut 

Then the functions uf satisfy the boundary value problems 

' uf'iy) = ut‘^{y) 

< ut i0) = 0 (6.18) 

, Uxiy/a±{x - te)) = y{x,e). 

Proof. This is an immediate consequence of the definition (6.17) and equation (6.15). □ 


Proof of Theorem 6.1. The differential equation in (6.18) is the same as in the Moranian case, 
so 

^tiy) = + y): 

where the period lattice has fundamental periods uf and The centering con¬ 

stant af is the larger zero of P^± (z) on (0, cuf). By the scaling laws for the P—functions, 

ut(y) = (ui -6 Afy) . (6.19) 

where = ujijujf. The boundary conditions at y = y/a±{x — te) in (6.18) and (6.19) imply 
that 

+ V^Xfix - te) = ( 6 . 20 ) 

oAx 
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By (6.13), \\mx^ooTl{xj£) = £, and furthermore x — fg ~ x, since = 0(1) as x — )• oo. 
Hence, —)• 0, by (5.9). Consequently, the right hand side of (6.20) goes to oo as x —)• oo. 
It follows that oi + ^/^A^x —)• cji, as x — oo. Therefore, 


lim A;^x = 

X—>oo 


uji — ai 


Cl 




( 6 . 21 ) 


Equations (6.17) and (6.19) imply that 

ut(y) = ut(v^y) = 6A^^p£i(ai + A^y/o^y). (6.22) 

Consequently, when 0 < y < x — 0 is fixed, Taylor expansion of Pci (z) around z = ai as 
in the Moranian case yields 

lim x^M^(y) = lim 6x^A^^P/:i(ai + y/o^A^y) 

X—^’OD X—^’CX) 

= 6(A^x)3y/EIPi^(ai)y+ 0(A^^) 

X—^’OD ^ 

(6£i) 6cfP^^(ai)y _ Ciy 

a± 

Similarly, if y = sx for some fixed 0 < s < 1, then by the continuity of Pjc^ (z) on (ui, tci) 
and (6.21), 

lim x^u^(sx) = lim 6(xA^)^P£^(ai + s-y/a±A^x) 

X—>-CxD X—^OO ^ 

(621) 6cfP£^(ai + sci) _ 02 (g) 
a± a± 

Here Ci and 02 ( 5 ) are as in Theorem 5.6. 

Finally, by Corollary 6.4, 

for each y fixed , —< lim x^Ux{y) < and 

a_|_ x->-oo a_ 

for each 0 < s < 1 fixed , x^Ua;(gx) < 

a_|_ x->oo a_ 

Letting <5 —)■ 0, so that a± we obtain (6.1) and (6.2). □ 


7 The Number of Killed Particles 


In this section we discuss the distribution of the number N = N^oi particles killed during 
the course of a critical branching Brownian motion with killing at 0 initiated by a single 
particle at position y > 0. Our primary interest is in the tail of the distribution, that is, in 
the large-A; behavior of the probabilities P^{Ai > fc}. 


Theorem 7.1. If the offspring distribution has mean 1, positive variance cr^, 
ment then 


m 

kF{N > fc) ~ 2Cjy\/rn where O 7 

k=l 



and finite third mo- 


(7.1) 
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The proof, which will use a form of Karamata's Tauberian theorem, will be given in 
section 7.4. 

If we knew that the sequence {N > k) were monotone then we could conclude from 
(7.1) that > /c) ~ C'^yjk^l’^. However, it seems unlikely that monotonicity of the se¬ 

quence k^y{N > k) holds in general. Thus, to obtain sharp asymptotic results about the 
individual probabilities P^(Ai = k), we will impose more restrictive hypotheses on the off¬ 
spring distribution that will allow us to avoid the use of Karamata's theorem. In its place, 
we will use a result of Flajolet and Odlyzko [FO90] that allows one to extract information 
about the asymptotic behavior of the coefficients of a power series from information about 
its behavior on the circle F of convergence. Our hypotheses are most conveniently formu¬ 
lated in terms of the functions 

h{s) = 2[T'(1 — s) — (1 — s)] and k{s) = f h{s') ds', (7.2) 

Jo 

where T' is the probability generating function of the offspring distribution. Since the 
power series for a probability generating function has radius of convergence 1, the func¬ 
tion h extends to an analytic function h{z) in the disk of radius 1 centered at z = 1 , as 
does its integral k. If the offspring distribution has finite support then the functions h and 
K are polynomials, and consequently are well-defined and analytic in the entire plane C. 
Observe that k has a zero of degree 3 at s = 0, since h{s) = + 0(|sp). 

Theorem 7.2. Assume that the offspring distribution has mean 1, positive variance , and that 
the function h{z) extends analytically to a disk of radius 2 +e centered at 0. If the indefinite integral 
k{s) has no zeros in the punctured disk 0 < |s| < 2, then for each y > 0, 

P^(A^ > A:) ~ where C 7 = (7.3) 

and 

]py(N = k) ^ ^ where Cg = (7.4) 

2v^ 

For the double-or-nothing (Moranian) offspring distribution, the functions h{z) and 
k(z) are given by h{z) = z^ and k,{z) = 2^/3, and so the conclusions of Theorem 7.2 
hold. (In this case, we will exhibit an explicit closed-form representation of the distribu¬ 
tion, in Theorem 7.9 below.) Consequently, by Rouche's theorem, the hypotheses of The¬ 
orem 7.2 hold for all finitely-supported offspring distributions in a neighborhood of the 
double-or-nothing distribution, in the following sense: for any integer m > 3 there exists 
am £ (0,1/2) such that the hypotheses of Theorem 7.2 hold for any probability distribution 
{qn}o<n<m such that min(g'o, ^ 2 ) > Om- 

Theorem 7.2 should be compared with recent results of Maillard [M13] and Berestycki 
et al. [BBHM15], which give sharp tail probability estimates for the number of killed par¬ 
ticles in the somewhat different context of supercritical branching Brownian motion with 
particle drift. Both [M13] and [BBHM15] also use the Flajolet-Odlyzko theorem, and so 
must also contend with the issue of analytic continuation of the generating function. In 
[BBHM15], the reproduction mechanism is simple binary fission, and so there is no need 
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to impose additional conditions. In [M13], the offspring distribution is arbitrary, but must 
have exponentially decaying tails and mean greater than 1. In all three cases, it would be 
of interest to determine optimal hypotheses on the offspring distribution. 

The proof of Theorem 7.2 will be given in section 7.5. In sections 7.1, 7.2, and 7.4, we 
shall assume only that the hypotheses of Theorem 7.1 are in force; in section 7.3 we shall 
assume that the offspring distribution is the double-or-nothmg distribution; and in sec¬ 
tion 7.5 we shall assume that the offspring distribution satisfies the hypotheses of Theorem 

7.2. 


7.1 Expected number of killed particles 
Proposition 7.3. 

= 1. (7.5) 


The proof will use the following a priori bound on the expectation. 

Lemma 7.4. 

Ey[N] < 1 

Proof. For each time t, let Nt denote the total number of particles frozen at 0 by time t, and 
Zt the total number of particles in the enveloping branching Brownian motion with no 
particle freezing (see the construction in section 2). The counting process dearly 

increasing in t, and 

lim Nt = N. 

Consequently, by the monotone convergence theorem, it will suffice to show that [Nt] < 
1 for any t > 0. 

Let Zt be the total number of particles at time t in the branching Brownian motion with 
freezing (including those particles frozen at 0), and Wt the initial particle's location at time 
i > 0. Define T to be the time of the first reproduction event (recall that this is a unit 
exponential random variable independent of the branching Brownian motion) and tq the 
first time that a particle reaches the origin. Since the offspring distribution has mean 1, 

pyzT^r, = 1 . 

Now a particle that reaches zero will, in the enveloping branching Brownian motion, en¬ 
gender a critical descendant branching Brownian motion, and so by the strong Markov 
property, 

IE°[-Z't_(TATo)] = 1- 

This implies that Ey[Zt] = ¥P[Zt] = 1. Clearly Nt < Zt, so E^[Ni] < E^[Zt] = 1 for any 

f > 0. □ 


Proof of Proposition 7.3. Let 


z; = Zt-Nt 


(7.6) 
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be the number of particles of the branching Brownian motion alive at time t that are not 
frozen. Since E^[Zi] = 1 for any t > 0, it is enough to show that EZ^ —^ 0 as t —)• oo. 

Clearly, z; < Zt, because the particles counted in are contained in the set of particles 

counted by Zt. In fact, 

Zt 

Z*t = ^ 1 {particle i trajectory C (0, oo)} 

i=l 

Zt 

E^Zl = E"^ ^ Ijparticle i trajectory C (0, oo)}. 

i=l 

To evaluate the last expectation, we use the discrete Brownian snake representation of 
branching Brownian motion described in section 2. Recall that in this construction parti¬ 
cles are represented by vertices of a Galton-Watson tree, and their locations are obtained 
by running conditionally independent Wiener processes along the edges. Thus, for any 
particle i counted in Zt, the conditional distribution of the trajectory {Wl}s<t of particle 
i up to time t, given the realization of the skeletal branching process, is that of Brownian 
motion started at y. Hence, 

Zt 

E^Zt = E^ ^ Ijparticle i trajectory C (0, oo)} 

= EZtP^{Ws does not hit 0 by time t} 

= py{Ws does not hit 0 by time t} 

—> 0 as f —)• oo. 


□ 


7.2 Probability Generating Function of N 

Define 

OO 

(p(y, s) = = k}s'^ (7.7) 

A:=0 

to be the probability generating function of the random variable N under the probability 
measure P^. Because the sequence is multiplicative. Proposition 3.1 implies that for 
each complex number s in the disk |s| < 1 the function y i->- ip{y, s) is and satisfies the 
differential equation 


dyy(p{y, s) = ip{y, s) - T'((p(y, s)) for y > 0. 

Furthermore, ip satisfies the boundary conditions 

p{0, s) = s and 
p{oo, s) = 1. 


(7.8) 


(7.9) 
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Because we are interested in the tail of the distribution, we will find it useful to re¬ 
formulate the boundary value problem for ip as an equivalent problem for the generating 
function 

CX) 

H{y, = (7-10) 

k=l 

Proposition 7.5. For each s such that |s| < 1, the function H{y, s) satisfies the differential equa¬ 
tion 

dyyH{y, s) = 

where h{z) = 2fF{l — z) — {l—z)]. In addition, H{y, s) satisfies the following boundary conditions: 

HiO,s) = s, (7.12) 

H{oo,s) = 0, and (7.13) 

\imH(y,s) = l. (7.14) 

Proof. The generating functions H and ip are related by 

OO 

ip{y,s) = J 2 ^y{N = k)s^ 

k=0 

OO OO 

= ^P^(Ar > A;)s^ - ^P^(iV > A: +l)s^ (7-15) 

k=0 k=0 

= 1 + H{y,s) (- 7 -) • 

Thus, the differential equation (7.11) follows directly from that for ip, as do boundary con¬ 
ditions (7.12) and (7.13). Finally, the additional boundary condition (7.14) follows from 
the hypothesis that the branching Brownian motion is critical, as this makes E^N = 1, by 
Proposition 7.3, and 

OO 

limH(y,s) = Vp 2 ^(A^ > k) =E 2 ^[iV]. 
k=l 

□ 


7.3 The Moranian Case 

Consider now the Moranian case, where the number of offspring is either 0 or 2, each 
with probability In this case the function h in the differential equation (7.11) reduces to 
h{z) = z"^, and so (7.11) becomes 

dyyH{y,s) = ^^H\y,s). (7.16) 

The boundary conditions (7.12) and (7.14) uniquely determine the solution, which can be 
written explicitly as 

H{y, s) = s ^-^yVl - s -6 1^ (7.17) 
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Remark 7.6. The function 


H{y, s) = s (- s - 1 


-2 


also satisfies the boundary value problem, but since it has a pole at y = y^6/(l — s), it 
cannot be a probability generating function for all y > 0. 


Similarly, let 

u{y,s) = 1 - ip{y, 1 - s), 

then the differential equation (7.8) and the boundary conditions (7.9) become 

dyyu{y,s) = u^{y,s) 
u{0, s) = s 
u{oo, s) = 0, 

which has the solution 

SO by (7.18), the generating function ip{y, s) is 


(7.18) 


(7.19) 


V’iy,s) = 1 - u{y,l - s) = 1- 


6 {l-s) 


(yVl - s + \/6)2 


(7.20) 


The equation (7.20) completely determines the distribution of N under P^. In Theo¬ 
rem 7.9 below, we will use (7.20) to provide explicit formulas for the probabilities P^ {N > 
k). But before doing so, we will show that (7.20) leads to the asymptotic formulas (7.3) and 
(7.4). 


Theorem 7.7. In the Moranian case, 

lim A:^P^(A^ > k) = 

and 


fc—>-oo 




■= c^y, 


lim kWy{N = k) = := Cey. 

k^oo 2v67r 


(7.21) 

(7.22) 


The proof will rely on the following theorem of Flajolet and Odlyzko [FO90]. 

Theorem 7.8 (Corollary 2, [FO90]). Assume that the power series A(z) = defines 

an analytic function in\z\ < 1 that has an analytic continuation to a Pacman domain 

Dob,s ■= {\z\ < 1 5} n {\arg{z - 1)| > /3} 

for some <5 > 0 and 0 < (3 < 7r/2. If 

A{z) ^ K{1 — z)'^ as z ^ 1 in Da, 5 , (7.23) 

then as n ^ oo, 

^ (7.24) 


F(-a) 


provided a ^ {0,1,2, • • • }. 
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Proof of Theorem 7.7. The functions H and (p given by equations (7.17) and (7.20) have al¬ 
gebraic singularities at s = 1, but for each y have unique analytic continuations to the slit 
plane C\{l<s<oo}. Expansion around s = 1 yields 


and 


H{y, s) = s/ - s + 1 

Vfe.s) = 1 - (1 - »)/ il + - 0)5 + ^(1 - ») 

= S + ;^{1-0)I+0(|1-S|"). 


(7.25) 


(7.26) 


Thus, the hypotheses of Theorem 7.8 are satisfied, and so the relations (7.21) and (7.22) 
follow. □ 


Because iE is a simple algebraic function of the argument s, its power series coefficients 
can be determined exactly. These provide an explicit formula for the distribution of N in 
the Moranian case. 


Theorem 7.9. In the Moranian case, 


where 


and 


^y{N>k) = Ak + Bk, 


„2k-2 


Ak = k{k + (?) 


Bk = 2c 


E 


i+j=k,i,j=0,l,- 


(c2 _ l)fc+l’ 

(2i - 3)!! 

¥i\ {? - 1 )J +1 ’ 


c = 




(7.27) 

(7.28) 

(7.29) 

(7.30) 


Proof The equation (7.17) for the generating function H{y, s) can be rewritten as 

H{y, s) = s(l cVl - s)“^- 
Expanding around s = 0 yields 

1 _ [1 - 


[l -|- C\/l “ [l -I- cVl — S ^ [l — C\/l “ s] ^ 

_ l-2cyT^+?{l-s) 

1 2 cyT^ 


+ c^ 


(l-s) 


[(1 _ c 2 ) + [(1 - c 2 ) + [(1 - C ^) + 

= I + II + III. 
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For I, we have 


/ = 


1 


(1 - c2) + C^S 

1 f 


1 


(1-C2)2 ^1-^ 
1 

(1-C2)2 
1 


(1-C2)2 

oo 

Y,ik + 1) 

k =0 





E 

j =0 


.2 _ 1 


„2fc 

^ h 

rs . 


(c2 - 1 )^+ 2 ' 


Now let 

then for III, we have 

III = 0^(1 -s)-I 

/ OO 

E 


ak = {k + 1) 


„ 2 k 


(c2 - 1)^+2’ 


oo oo 

2 I _fc _ k +1 


= c“ I 7 ^OfeS- - 


^A:=0 fe=0 


= c 


= c 


'^aks’^ - 
/c =0 fc=l 

oo 

no + ^{ak — ak-i)s 


k=l 


(731) ^2 ( 1 


+ E 


(Pk -2 


(c2 _ 1)2 ' ^ (c2 _ l)fc+l 


^ 2 k 


E 

k =0 


( c 2 - 1)^+1 [ c 2 - 1 


-ik + l)-k 


c2 — 1 




{k + l)-k 


Next, set 


„2fe 


Pk — 


{k + l)-k 


(c2 - 1)^+1 [c2 - 1 
then we find that Ak+i = ctk + 13k by direct computation 
It remains to show that 


II = J2Bk+ii 


k =0 


(7.31) 


(7.32) 
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By Newton's binomial formula. 


Consequently, 




i=0 


II = -2^ ^ 


1/2 


-S 


vi=0 

' °° 1 Mi-1/ 


E 

vi=o 




Oil S' 


IE' 


ii=0 


M=0 


24E 


ii=0 


(2z-3)!! 
2^i! ' 


E 


O/S-' 


2-E E 

/c=0 \ 2+_7=fc,^j=0,l- 


yj=0 

(2i-3)!! 


2*i! ^ 


ai s 


( 7 . 31 ) 


^ ^ Bk+li 


k=0 


(7.33) 


□ 


7.4 Proof of Theorem 7.1 

In the general case, where the offspring distribution is assumed only to have mean 1, vari¬ 
ance > 0, and finite third moment, the behavior of the generating functions (p{y, s) and 
H{y, s) as s —1— can be deduced from that in the Moranian case by comparison argu¬ 
ments, as in section 6. This derivation exploits the fact that for s near 1 the differential 
equation for H looks like that for the Moranian case, for which we have exact solutions. 
The result, the details of whose proof we defer to section 7.6, is as follows. 

Lemma 7.10. For each y > 0, the generating functions (p{y, s) and H{y, s) satisfy 

ip{y, s) — s ~ -^(1 — and (7.34) 

V6 

H{y,s)-l^^{l-s)^/^ (7.35) 


flS S t 1- 

Because the generating functions ip{y, s) and H{y, s) are defined by power series with 
nonnegative coefficients, the singular behavior of their derivatives can be deduced from 
Lemma 7.10, by the following elementary fact. 
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Lemma 7.11. Let A : [0,1] —)• M+ be an absolutely continuous, nonnegative, increasing function 
whose derivative A' is non-decreasing on (0,1). If for some constants C > 0 and a G (0,1), 

^(1) - A(s) ~ C(1-s)" fls stl, (7.36) 

then 

A'(s) ~ Ca(l - as s 11- (7.37) 

Proof Since A is absolutely continuous, 

PSl 

A(si) — A(so) = / Lift) dt for all 0 < sq < si < 1- 
J So 

Suppose that for some (5 > 0 there were a sequence Sn —)• 1— along which A'{sn) < Ca{l — 
(5)(1 — Since A' is non-decreasing, it would then follow that < C'q!(1 — (5)(1 — 

for all s < Sn, and so for any e > 0, 

- A{sn{l - e)) < Ca{l - 5)(1 - Sn)°‘~^e. 

But this would lead to a contradiction of the hypothesis (7.36) provided e is sufficiently 
small relative to 6. A similar argument shows that it is impossible for A'{sn) > Ca{l + 
(5)(1 — s„)““^ along a sequence ^ 1—. □ 

Corollary 7.12. For each y > 0, as s ^ 1—, 

1 OO 

-H{y, s) = ^ kFy{N > k)s^ ~ ^(1 - (7-38) 

Theorem 7.1 follows directly from Corollary 7.12 and Karamata's Tauberian theorem 
(cf. [BGT87], Corollary 1.7.3), which we now recall. 

Theorem 7.13. Let A(z) = be a power series with nonnegative coefficients an and radius 

of convergence 1. If, for some constants C,(3 > 0, 

^(s) ~ C/(l-s)^ flsstl, (7.39) 

then as n ^ oo. 

n 

J2ak-^CnP/r{l + f3). (7.40) 

k=l 

□ 


7.5 Proof of Theorem 7.2 

In this section we assume that the offspring distribution satisfies the conditions enumer¬ 
ated in Theorem 7.2, in particular, that the function k{z) defined by (7.2) has no zeros z 
such that 0 < \z\ < 2. We will once again make use of the Flajolet-Odlyzko theorem 
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(Theorem 7.8), which requires (a) that the function defined by the power series in question 
should vary regularly as functions of s near the singularity s = 1, and (b) that this function 
has an analytic continuation to a Pacman domain. Lemma 7.10 implies that <p(y, s) and 
H{y, s) vary regularly as s —)• 1— along the real axis from below. The following lemma 
ensures that they have analytic continuations to a Pacman domain, and that the regular 
variation persists in this region. 

Lemma 7.14. If the offspring distribution satisfies the hypotheses of Theorem 7.2, then for each y > 
0 the generating functions H{y, s) and y>{y, s) have analytic continuations H{y, z) and (f{y, z) to 
a slit disk {\z\ < 1 + (5} \ {1 < z < 1 + 5}, and the relations (7.34) and (7.35) hold as s ^ 1 in 
the slit domain. 


Proof. In view of the relation (7.15), to prove that the function H{y,s) has an analytic contin¬ 
uation it suffices to show that the probability generating function y>{y,s) can be analytically 
continued, or alternatively that the function 

u{y,z) := 1 - <p(y, I- z) 


has an analytic continuation to a slit disk {|1 — zj < 1-|-5}\{1<1 — 2 ;<l-|-(i}. Recall 
(Proposition 3.1) that for all real z in fhe inferval 11 — ^| < 1 the function u{y, z) satisfies the 
boundary value problem 

dyyu{y,z) = h{u{y,z)), 

u{0,z)=z, (7.41) 

u{oo, z) = 0. 

Integration of this differential equation, as in section 3 of [SF79], leads fo fhe equation 


/ 


dy 


{x,z) y/nffj) 


= x where K{y) := 2 / h{y')dy'. 


(7.42) 


(The upper limit of integration is z because m(0, z') = z.) 

It is easily checked that in any region of the z—plane where the equation (7.42) has 
a solution u{x, z), the solution will satisfy fhe boundary value problem (7.41). Thus, to 
prove the first assertion of the lemma it will suffice to show that the integral equation 
(7.42) implicitly defines u{x,z) as an analytic function of z for 2 ; in a slit disk. 

Define a function of two complex variables z, w by 


G{w, z) 



dy 


This function is analytic in z and w in any domain D C such that there exists a simply 
connected domain D C C in which 1 / ^/K is analytic and such that for any z,w ^ C fhere is 
a path from z to w m D. Moreover, 


dz 


l/^/^ and 


dw 


l/y/K{w). 
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Therefore, by the complex implicit function theorem, the equation G{w{z), z) = x defines 
w{z) as an analytic function of z in a neighborhood of any solution G{wo, zq) = x where 
+ 0 . 

The function <p(x, z) is analytic in the unit disk and, since it is a probability generating 
function, satisfies |(p| < 1. Consequently, the function u{x, z) = 1 — ip{x, 1 — 2 ) is analytic 
in the disk Di := {\1 — z\ < 1} and satisfies tt(x, z) G Di for all z G Di. By hypothesis, the 
functions h{z) and k{z) have analytic continuations to a disk of radius 2 + 5 > 2 centered 
at 0, and the only zero of k{z) in this disk is at 2 : = 0. Therefore, the functional equation 

G(u(x, z), z) = X (7A3) 


holds for all z £ Di. 

We claim that for any point ^ G dDi except ^ = 0, the function u{x,z) converges as 
2 ; —^ to a value u{x, such that u{x,^) G Di. This is clearly equivalent to the assertion 
that for any point / 1 of the unit circle, the function (p{x,z) converges as 2 : ^ to a 
value (p{x, of absolute value less than 1. To see that this is so, recall that (p{x, z) is the 
probability generating function of the random variable N under P^. It is easily seen (for 
instance, using the discrete Brownian snake construction) that for any x > 0, 

P^{N = A:} > 0 for every A: = 0,1, 2 ,... . 

But this implies that is less than 1 for every 0 G [—vr, tt] \ {0}. 

It now follows that the functional equation (7.43) extends by continuity to all 2 ; G dDi 
except 2 = 0, and that at any such boundary point, u{x, z) / 0. Hence, the function 
1/ ^/K (w) is analytic and nonzero in a neighborhood ofw = u{x, z). This implies that u{x, z) 
has an analytic continuation to a neighborhood of every 2 G dDi except the singular point 

2 = 0. 


Next, we must prove that u{x, 2 ) has an analytic continuation to a slit neighborhood of 
2 = 0. Since (p(x, 1) = 1, the function u{x, 2 ) converges to 0 as 2 0 in the disk {| 1 — 2| < 1}; 
thus, the function k{u) approaches 0. Our assumptions on the offspring distribution imply 
that 

CXD 

h{z) = cr^ 2 ^ + ^ hjZ^ and hence 
i=3 

2 °° 

Kiz) = -a^z^ + ^ A:j2^; 

i=4 


consequently. 


^==iT2-3/2(l + ii(^)) 
^/k{z) 


where K = \j3j2^ and R{z) is an analytic function in some neighborhood of 2 = 0 and 
satisfies i?(0) = 0. Now the integral equation (7.42) and the implicit function theorem 
imply that, for fixed x > 0, the function u = u(x, 2 ) satisfies the differential equation 


du _ a/ k{u) (1 + R{z)) 

dz ^k(z) z^/"^ (1 + R{u)) 
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for all z in a domain {|1 — 2 :| < l}n{| 2 :| <5}. Using the analyticity of R, we conclude that 

^ bkU^ = ^ + c, 

where C is a constant of integration. Squaring both sides exhibits u as a meromorphic 
function of -y/z. This shows that u has an analytic continuation to a slit disk {| z| < 5} \ 
{—6, 0], and so it follows, by the relations u{x, z) = 1 — if(x, 1 — z) and equation (7.15), that 
(f(x, z) and H{x, z) admit analytic continuations to a slit disk {|1 — 2 ;| < 5'} \ [1,1 + 5']. 

Finally, we must prove that the relations (7.34) and (7.35) hold as z —)• 1 in the extended 
domain of definition of fhe functions <p(x, z) and ff(x, z). But the analytic continuation 
argument above shows that, in a slit disk centered at z = 1, the functions (p(x, z) — z and 
H{x,z) — 1 are meromorphic functions of y/1 — z, and hence have Puiseux expansions in 
powers of (1 — zY^^. Since (7.34) and (7.35) hold as s t 1, it follows that they persist in the 
slit disk. 

□ 

Proof of Theorem 7.2. Lemma 7.14 implies that the generating functions (p(y, s) and H{y, s) 
meet the requirements of fhe Flajolet-Odlyzko theorem (Corollary 7.8). Therefore, Theo¬ 
rem 7.2 follows from relation (7.24) (since r(—i) = —2y/7r and r(— |) = |y/vr)- 

7.6 Proof of Lemma 7.10 

The strategy is similar to that of section 6. As s —)• 1, 

«(y, 1 - s) 0 

- — -H{y,s) 0, 

s 

and so for s near 1 the differential equations (7.41) and (7.11) for u{y, 1 — s) and H{y, s) 
have forcing terms that are nearly quadratic. Taylor expansion of h shows that for any 
(5 > 0 there exists e > 0 such that for 1 — e < s < 1, 


a-u‘^{y, 1 — s) <h{u{y, I — s)) < a+u^{y, 1 — s) and 



1 — s 

^ /l — s \ 

1 — s 

a_ 

- B{y,s) 

s 

<hi ^ B{y,s)j < a+ 

- B{y,s) 

s 


where a± are defined in (6.14). Let u±{y, 1 — s) and H±{y, s) satisfy the boundary value 
problems 


' dyyU±(y, 1 -s) = a±u±(y, 1 
< u±(0, s) = s 

u±(oo,s) = 0, 


s) 


f Q 

dyyB±{y, s) = a±—^B\{y, s) 

' iF±(0,s) = s 

, iT±(cx),s) = 0, 
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and set ^±{y, s) = 1 — u±{y, 1 — s). By the same argument as in Corollary 6.4, 


u+{yA-s) <u{y,l- 
(p-{y,s) <(p{y,s) <_ 
H+{y,s) <H{y,s) 

Define re-scaled versions 

u±{y,l - s) = u±{ 
H±{y,s) = H± 

these satisfy the boundary value problems 


_ _ g _ 

dyyH±{y,s) = 

H±{0,s) = s 
H±{oo, s) = 0, 

which are the same as in the Moranian case. Hence, (p±{y, s) = l — u±{y,l — s) and H±{y, s) 
have the same asymptotics as (7.26) and (7.25): 

f±{y,e) “ * + - s)5 +0 (|1 - sp) 

as s —)• 1. Hence, applying (7.45) yields 

^±{y,s) = s + {1 - s)^ -6 0(^|l-sp) 

H±{y,s) = 1 - +0(|1 - s|). 


' dyyU±{y, 1 - s) = n±(y, 1 - s) 

< w±(0) s) = s < 

u±(oo,s) = 0, 


s) < u-{y,l - s) 
^+{y:S) 

< H-{y,s). 


,1- s 


(*'•) 


(7.44) 


(7.45) 


Finally, as 5 —)■ 0 we have a± , and so (7.44) implies 

‘lay 3 

(p(y,s) - s ~ -^(1 - s )2 

Hte.s)-1-^(1-s)i 


as s —)• 1. 


□ 
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